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We study a twisted Hubbard tube modeling the [CrAsJoo structure of quasi-one-dimensional su¬ 
perconductors A. 2 Cr 3 As 3 (A = K, Rb, Cs). The molecular-orbital bands emerging from the quasi¬ 
degenerate atomic orbitals are exactly solved. An effective Hamiltonian is derived for a region where 
three partially filled bands intersect the Fermi energy. The deduced local interactions among these 
active bands show a significant reduction compared to the original atomic interactions. The re¬ 
sulting three-channel Luttinger liquid shows various interaction-induced instabilities including two 
kinds of spin-triplet superconducting instabilities due to gapless spin excitations, with one of them 
being superseded by the spin-density-wave phase in the intermediate Hund’s coupling regime. The 
implications of these results for the alkali chromium arsenides are discussed. 

PACS numbers: 71.10.Pm; 72.15.Nj; 74.20.Mn; 74.70.-b 


Introduction. —Recently, the alkali chromium arsenides 
A 2 Cr 3 As 3 (A = K,Rb,Cs) have been found as a new fam¬ 
ily of inorganic quasi-one-dimensional (Q1D) supercon¬ 
ductors with strong electron correlations^ 3 The basic 
building block of these compounds is the [CrAs]6 clus¬ 
ter consisting of two conjugated triangular complexes 
[CrAs] 3 as shown schematically in Fig Oja). They are 
aligned along the c axis forming a [CrAsJoo tube, and 
intercalated by A + cations forming a hexagonal lattice. 
The density functional theory (DFT) calculations^ pre¬ 
dict a three-dimensional (3D) Fermi surface (FS) sheet 
(7 band) and two Q1D FS sheets (a and /3 bands), essen¬ 
tially due to the Cr 3d electrons. The NMR experiment^ 
has revealed a power law behavior of the spin-lattice re¬ 
laxation rate, manifesting the Luttinger liquid feature 
above T c . The penetration depth measurement^ has ev¬ 
idenced a line nodal feature in the pairing state below 
T c . 

Because of the existing 3D 7 band, whether the super¬ 
conductivity is solely originated from the Q1D structure 
of A 2 Cr 3 As 3 is uncertain. In fact the nearly isotropic 3D 
bulk Cr As compound shows superconductivity with T c ~ 
2.2 K under pressure of ~ 0.7 GPa££ Interestingly, Zhou 
et al. pointed out that an /-wave pairing state could arise 
from the 3D band with a node line while a fully gapped 
p-wave pairing state could dominate at the Q1D bandJ^ 
Such triplet superconductivity, with some variations in 
spatial symmetryii, could be driven by ferromagnetic 
fluctuations within the sublattice of Cr atoms^ii. 

In order to understand the formation of the low energy 
bands, it is particularly important to understand the elec¬ 
tronic property of a single fundamental [CrAsJoo tube. In 
this Letter, we model this system by a twisted Hubbard 
lattice composed of triangular complexes coupled along 
the c axis with the glide reflection symmetry as shown 
in FigQJb)^ In each unit cell there are six Cr atoms, 
each with five 3d atomic orbitals (AOs). The influence 



FIG. 1: (Color online) (a) A CrAs cluster in the ab plane. 
The solid (dotted) circles connected by the solid (dotted) lines 
represent the Cr atoms in the first (second) triangle in a unit 
cell. The isolated outer solid (dotted) circles represent the As 
atoms in the corresponding planes, (b) A Q1D CrAs tube. 
The blue (green) filled circles represent the Cr atoms in each 
triangles. The As atoms are not shown. 


of the As 4 p orbitals can be effectively accounted for the 
indirect hopping of Cr-3d electrons. So the model in¬ 
volves thirty energy bands in total. In the realistic case, 
fortunately, only three partially filled bands are active 
in the low energy regime. We will explicitly show how 
these bands come from the molecular orbitals (MOs) of 
[CrAs]gMi. Our purpose is then to understand their co¬ 
operative low temperature physics within the Luttinger 
liquid approach. The proposed effective model is of in¬ 
terest in its own right as we shall explore in the following. 

Model Hamiltonian .—The Hubbard model for a single 
[CrAsJoo tube is expressed as H = Ho + H [ nt , where, 
Ho represents the noninteracting part consisting of the 
tight-binding kinetic energy and the crystalline electric 
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field (CEF) splitting, 


Ho = - t^ r Jdl^d^+ Y J E rm n^.{l) 


rr'mm' a 


Here, d denotes the annihilation operator of Cr 3 d 
electrons at the site r with angular momentum m = 
0 , ± 1 ,± 2 , spin a =t,|- nffi and are the corre¬ 

sponding density and spin operators. The two twisted 
Cr triangles could have different E rrn = E^\ Em\ ac¬ 
countable for the possible occupation difference^, while 

rm rm/rn'cra' 

- J H Y, S r m) • Sr™ 0 (2) 
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represents the local interactions including the intraorbital 
Coulomb interaction [7, the interorbital Coulomb inter¬ 
action t/', the Hund’s coupling J#, and the pair-hopping 
Jp, respectively. 

There are four kinds of adjacent intraorbital hoppings 
trnm = tm (i = 1 — 4), corresponding to the nearest- 
neighbor (NN) sites in the first and second triangles, and 
those between the intracell and inter cell triangles, respec¬ 
tively, as illustrated in FigOJb). Because of the metallic 
bonding among Cr atoms, the direct orbital mixings are 
relatively small, and the indirect hybridization is mainly 
bridged by the As 4 p orbitals. So it is legitimate to con¬ 
sider a simpler situation for the adjacent interorbital hop- 

P in g ; tmm> = VtmS | m |,|m'| for TO ^ TO 7 , with \r]\ < 1. In 
this situation, the atomic orbitals are quasidegenerate as 
the nonvanishing mix terms are isotropic in spaced. Fi¬ 
nally, we include the next NN intraorbital hopping t$ 
along the tube direction. 

Molecular-orbital bands. —Denoting each site by 
r = (n, a, £), with a = 1,2,3 being the location in the 
first (£ = 1 ) or second (£ = 2 ) triangles in the nth 

unit cell, it is convenient to introduce a base d^ = 

\ T 

J for 

the atomic m-orbital in the nth unit cell (the spin index 
a is implied). For m = 0, this base accommodates a 
representation for the C 3 rotational symmetry, leading 
to six MOs corresponding to E,E',A, and A' states, 
respectively^. For mm E 1 or ±2, we need to introduce 
a set of new base = ^[d^ A Thus for 

a single [CrAs ]6 cluster, we have thirty MOs defined by 

Cn ^ = (^0 Qo)dn^ for r = 0 (denoting dn^ = dn°^), ±1 
and ± 2 , respectively, with (uj = e lcp , cp = 2tt/3) 
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Note that the eigenstates of R with eigenvalues Ai = 2 
and A 2 = A 3 = — 1 constitute of representations A and E 
(or A! and E’\ respectively. 

When the triangles are coupled along the c axis via the 
intercell hopping t ^ m ,, we can extend Qo to the momen¬ 
tum k- resolved matrix Q^a\k) so that the Bloch form 

M t ) — (A r ) A r ) A t ) A t ) A t ) A t ) \ t 

^k 1 , 1 ) ’ (/e,l, 2 ) ’ C (/c, 2 ,l) 5 c (/c, 2 , 2 ) ’ C (/c, 3 ,l) ’ C (/c, 3 , 2 ) J 

is still a natural base diagonalizing i7 0 , leading to thirty 
MO bands labeled by the eigenenergies £^^(k). Here 
the subscript £ = 1,2 corresponds to the antibonding or 
bonding bands, respectively, due to the twisted structure. 
The explicit expressions of Q^a\k) and £^^(k), which 
also depend on the orbitals r(= 0 , ± 1 , ± 2 ) and C 3 eigen¬ 
values A a (a = 1 , 2 , 3), are presented in the Supplemental 
Material (SM)i£. A set of subscripts (a,£) determines 
the symmetry property of the corresponding MO bands. 

We fitted the DFT band structure along the tube di¬ 
rection using the obtained MO bands within r = 0 , ± 2 , 
while the bands with r = ±1 are fairly away from the 
Fermi energy as revealed in the DFT calculations^. The 
three partially filled DFT bands, i.e., the 3D 7 band 
characterized mainly by the d z 2 orbital (m = 0 ), the 
Q1D a and /3 bands characterized mainly by the d xy and 
d x 2_ y 2 orbitals ( \m\ = 2 ), are all holelike near the T 
point (fc = 0) and electronlike near the A point (k = 7 r). 
Therefore, the 7 band corresponds to the singlet MO 
band labeled by (r = 0,a = 1,£ = 1). The a and j3 
bands, which are degenerate along the whole T —>> A di¬ 
rection, correspond to the doublet MO bands labeled by 
(r = —2, a = 2,£ = 2) and (r = —2, a = 3, £ = 2), re¬ 
spectively. The best fitting using £(i\y £( 22 ) = ^( 32 ) 
is shown in Fig|2|^. Here, the tight-binding parameters 
are not uniquely determined because the number of these 
parameters exceeds eight necessary coefficients in the fit¬ 
ting. On the other hand, the precise values of the fitting 
parameters are not important in the present study. As 
we shall find later, only symmetry property of the MO 
bands and local interactions between them play a crucial 
role in the Luttinger liquid approach. 

For simplicity, from now on, we shall use the band sub¬ 
script v (= 1,2,3) to account for the three MO bands in¬ 
tersecting the Fermi energy. These active MO bands are 
associated with phases p v = 0, 27r/3, — 27r/3, or chiralities 
d v = 0, 1, and —1, respectively. In the full ID Brillouin 
zone, there are three pairs of Fermi points (kp v% — kp u ), 
satisfying 0 < kp 1 < kp 2 = kp 3 < 7 r, as schematically 
shown in the inset of Figj2j By integrating out all inac¬ 
tive bands, we obtain the effective theory describing the 
low energy property of the active bands: 


Heff — ^ ^ £v{k')'dkv<j T 


(n) 
int ' 
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( 3 ) Where, n^ va = c k ua c ^o- is the density operator of elec¬ 
trons in the zAh MO band, the residual short-range 
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FIG. 2: (Color online) Fitting the band structure: The upper 
band is twofold degenerate. The lower inset is the schematic 
picture for the three partially filled bands with three pairs of 
Fermi points. 


interactions in the MOs in the nth unit cell, given by 
H int = + T U vv 'h ua (n)h v ’ a '(n) 

V v^v'tjcr' 


- y. J vv >$>v(n) ■ s„/(n) (5) 

v^v' 

+ Ji23 4 t (n)cJ i (n)c 2 4 .(n)c 3 t (n) + (2 3) + H.c. 


In this expression, the electron annihilation operator in 
the n-th unit cell is defined by the Fourier transformation 
cw(n) = Yl,k elkncac kvcj (with c 0 the lattice spacing 
taken as unit). Only those terms preserving the neu¬ 
trality condition = 0 could survive. The ma¬ 

trix Qo is used in deducing Eq.(5) as the short-range 
interactions are mainly due to the slowly varying part, 
leading to IJ\ = [7/6, U2 = U3 = (U + U' + Jh + J p )/ 12 ; 
U 12 = U 13 = U'/12-J h /24; U 23 = (U+U'+J H +J P )/ 48; 
J12 = J13 = J23 = (?7 + U' + Jh + « 7 P )/ 12 ; and 

-[123 = Jp/6. The influence of inactive bands is mainly 
accounted to the renormalized tight-binding parameters. 

The Luttinger liquid in the weak-coupling regime .— 
We now take the continuous limit, linearize the active 
bands near the Fermi points, and decompose the elec¬ 
tron operator into right and left moving components like 
Cva(z) ~ e-^^-^L^z) + Here, 

z = nco is the spatial coordinate along the tube direc¬ 
tion, R V(T and L vc r represent the right and left moving 
fermions describing the low energy excitations near the 
Fermi points (kp u , —kp u ) with linear dispersion ±vp u k. 
The long-wavelength, low-energy effective Hamiltonian 
(density) is given by H e ff = Ho + 'Hint, where Ho = 
J2u,A iv F^)l R t<T d zRua - Ll a d z L„ a ] is the kinetic part, 
and Hint includes various residual interactions which are 
usually expressed in terms of the g-olog y 16 ! 17 . We shall 
assume the Fermi velocities vp u to be the same as this 
does not influence the nature of superconductivity we 
concern. The corresponding one-loop renormalization 


group (RG) equations resemble those for three-leg Hub¬ 
bard ladder s 18 ! 19 or a variant of carbon nanotubes 22 r—. 
The instabilities of these RG equations are classified rou¬ 
tinely: (i) the intraband instabilities as those developed 
in the single-channel Luttinger liquid^, and (ii) the inter¬ 
band instabilities as those developed in the two-channel 
band Luttinger liquid. Note that the three-band interac¬ 
tion in Eq.(5) does not lead to the peculiar three-band in¬ 
stability suggested in Ref^ as shown in the S M 15 i 24 . All 
these suggest the validity of the conventional bosoniza- 
tion approach based on spin-charge separation, where 
various ordering instabilities can be determined by Lut¬ 
tinger parameters. The new ingredients here are the pe¬ 
culiar symmetry surviving in the active MO bands and 
their dependence on local electron interactions. 

The right- and left-moving fields are then expressed in 
terms of the charge fields ( 0 ^ c , 0 ^ c ) and the spin fields 
(0i/, aj 0i/,s) (for each v = 1 , 2 , 3) by 

a (z) = ^ R,UC7 c i \/' K /‘2(0c.v+v0s.v-4>c.v-(T(l>s.v) 

V 2 ttc 0 

L v Jz) = FL,va P »yV/2(e(g) 

’ V27TC0 

The Klein factors Fh, vg and Fl^ g ensure the fermionic 
statistics between the right and left moving fermions. 
Next, in order to diagonalize the kinetic part, we need 
to introduce a set of new base 


r7,i 


^ 7,2 (^7,i07,l “1“ 07,2 T 07,3) , 


<7,3 = ~^J= (-07,2 + 07,3) : 
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where 7 = s, c, q 1 7 

2 Fl 2 L _ 2t/23 

_ 1 u r — _ 


= for i = 1 , 2 , 


Tr’ “ 2 ?F> are the 


a c = be = a s = b s = 

normalization constants. Similar relationships apply to 
the fields and for i = 1 , 2 , 3. Hence, we arrive at 
the following three-channel Tomonaga-Luttinger liquid 
Hamiltonian: 


*-/ 


dz Y 

i=l,2,3, r y=s,c 


|^-^-( V ^7,i) 2 + V,i(V0 7 ,i) 2 j (8) 


where, A= t 1 + 
and 3 = t 


& 7 - (- 1 )’ 


8 a 7 


6 7 


for i = 1 , 2 , 


a 7 ,3 = r 7 — 6 7 , t c = if+ |7 t a = Therefore, 

the Luttinger parameters are obtained explicitly by 


Kr 


K s ; = 


l+4^7-Sv / 8(2t7-5J ff )2 + C/2 


1 - 




_ u _ 

47 tvf 12ttvf 


for the channels i = 1 , 2 , respectively, and K C: 3 = K s ,3 = 
1 for the third channel i = 3. Here, we have adopted 
the conventional relations J p = Jh and U' = U — 2 J#, 
reflecting the rotational symmetry of the original AOs 2 £. 

Now since U > 0 and Jh > 0 , one can find that: (i) 
K c , 1 < 1 in the entire region and K S: 1 < 1 only when U < 
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FIG. 3: (Color online) Luttinger parameters for the channel-2: 
K s ,2 > 1 everywhere, and K c ,2 < 1 only in the intermediate 
regime between the lines Jh = 0.2 U and Jh = 0.6LL 


Jh ; (ii) K C 2 < 1 in the region 0.2 U < Jh < 0 .6U and 
K S 2 > 1 in the entire region. Specifically, in the physi¬ 
cally relevant regime, U > Jh, the spin excitations are al¬ 
ways gapless, so K s ^ could be fixed to the unit due to the 
spin-SU( 2 ) symmetry. Because K C: i < 1 , the channel- 
” 1 ” is in the spin-density-wave (SDW) phas o 16 ! 17 . The 
channel-” 3” involves the antibonding of the MO bands 
v — 2,3 as shown in Eq.([7J). In this channel both spin 
and charge excitations are critical. Because of the ab¬ 
sence of a spin gap, the dominating superconducting in¬ 
stability is the interband spin triplet pairin g 16 ! 17 , driven 
by the interband scattering between the two Q1D a and 
P bands. The intriguing case is the channel-” 2”, whose 
property depends on the ratio Jh/U. We plot in Fig. [3] 
the phase diagram determined by the Luttinger parame¬ 
ters in this channel. We find that K c 2 > 1 in the regimes 
separated by the orange-dotted and blue-dashed lines, re¬ 
spectively. In these two separated regimes, the dominat¬ 
ing instability is still the spin-triplet pairin g 16 ! 17 . But in 
the intermediate regime where K C: 2 < 1, the SDW insta¬ 
bility dominates. It should be noticed that in either case 
where the interband triplet superconducting instabilities 
dominate, the spin-singlet superconducting instability is 
the subdominating instabilit y 17 ! 26 . 

In order to see whether the above Luttinger liquid 
results are robust against deviations from the atomic 
orbital rotational symmetry, we have also considered a 
small deviation AU away from the rotational symmetry 
by assuming [/' = [/ ± AU — 2 Jh- As shown in the SM, 
the channel-3 is still in the critical phase and the role of 
AU is to modify the value of U in a simple manner so 
that the results remain unchanged^. 

Summary and discussions. —We have focused on the 
microscopic formation of the MO bands in a twisted Hub¬ 
bard tube capturing the Q1D nature of K^CrsAss, a new 
Q1D multiorbital superconducting molecular crystal with 
the moderate Coulomb interaction and Hund’s coupling. 
A three-channel Tomonaga-Luttinger liquid Hamiltonian 
describing the low energy physics of the three active MO 
bands ( the O', /?, and 7 bands) is then derived, showing 
possible unconventional triplet superconducting instabil¬ 
ities within a reasonable range of interaction parameters. 

The conclusions and implications of our study are com¬ 
pared with the previous studie s 10 ! 11 where a phenomeno¬ 


logical 3D Hubbard model for the three active MO bands 
was proposed based on the elegant symmetry argument^ 
and investigated by the random phase approximatio n 10 ! 11 
and the mean field treatment^. First, the twisted struc¬ 
ture of the Q1D [CrAsJoo tube showing the extended glide 
reflection symmetry in accordance with the C 3 group is 
explored in our approach so that the symmetry property 
of all the thirty MO bands (including the three active 
MO bands) could be identified. Second, the interactions 
among the three active MO bands are derived from the 
microscopic atomic Hubbard interactions, different from 
those proposed phenomenologically. Third, the three di¬ 
agonal channels in our Tomonaga-Luttinger Hamiltonian 
are superpositions of the original DFT bands, in contrast 
to the random phase approximation approach and the 
mean field treatment. 

We found two kinds of spin-triplet pairing instabili¬ 
ties emerging out from two of the three channels. One 
involves the Q1D a and ft bands, another involves all 
three bands. In the Luttinger liquid approach, triplet 
pairing instabilities are due to gapless spin excitations 
for U > Jh and K c > 1 in the corresponding chan¬ 
nels. The ferromagnetic correlation within the sublattice 
of Cr atoms^iii, though possible, is not a prerequisite of 
the triplet states. We also found an intermediate regime 
0.2 < Jh/U < 0.6 where the SDW phase emerges. Our 
solution is sensitive to the symmetry or regularity of the 
two conjugated Cr triangles, seemingly consistent with 
the recent hydrostatic and uniaxial pressure experimen¬ 
tal studj*2I. The exact mapping from the AOs to MO 
bands will also pave the way for further investigations on 
related effects such as the spin-orbit coupling within a 
microscopic framework. 

As the present study is limited to the Q 1 D case, the 
spatial symmetry of the superconducting pairing states is 
not specified. However, the actual 3D superconductivity 
can be perceived based on the Q1D physics because the 
identification of the three low energy MO bands is robust 
owing to the same symmetry argument. The local inter¬ 
actions among the MO bands are similar to those in the 
Q 1 D case. If the local atomic interactions are estimated 
as those in other Cr-based oxides, like SrCrO^, one has 
U ~ 2.7 ±0.5 eV, J H ~ 0.42 ± 0.1 eV, and J H /U ~ 0.16, 
then in the Q1D case the channel-2 is in the triplet phase 
but close to the SDW low boundary Jh/U = 0.2 shown 
in Figj3j The corresponding residual MO interactions U 
and J in the 3D case are significantly suppressed, but 
the ratio J/U enhanced-^, corresponding to the regime 
with small U but relatively large J/U in Ref^, where 
the triplet f y (^ x 2 -y 2 ) pairing state is favored. Of course, 
we have not considered the long-range Coulomb interac¬ 
tion and the electron-phonon coupling, the suppression 
of residual MO interactions should necessitate further in¬ 
vestigations on these influences. 

Finally, a more intriguing issue is the possible dimen¬ 
sional crossover from Q1D to 3D which could be tuned by 
either chemical substitution^ or physical pressure 27,29 . 
On one hand, one of the three active bands, correspond- 
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ing to v — 1, evolves with the intertube hopping and 
crossovers to the 3D 7 band which could lead to the line 
nodal feature. Meanwhile, the (<a, /3 ) bands could remain 
in Q1D because the intertube hopping among the AOs 
with m = ±2 is reasonably small. On the other hand, 
the Q1D superconducting instability can lead to a true 
long-range order when the intertube hopping is taken 
into account. Recall that the interband triplet pairing 
instability in the channel-2 is driven not only by the 7 
band, but also by the (<a, (3) bands. Consequently, the 
spin-triplet pairing instability in the channel-2 involves 


both the 3D and Q1D bands. As such a 3D pairing state 
could emerge from a normal state of an essentially Q1D 
Luttinger liquid characteristic, a scenario which is likely 
consistent with available experiments. It is desirable to 
investigate the related crossover behavior in this class of 
materials in the future. 
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By: Hanting Zhong, Xiao-Yong Feng, Hua Chen, and Jianhui Dai 


SUPPLEMENTAL MATERIAL 


This is the Supplemental Material (SM) for our paper titled ’’Formation of Molecular-Orbital Bands in a Twisted 
Hubbard Tube: Implications for Unconventional Superconductivity in K^CrsAss” i In this SM, we provide the detailed 
solution of the tight-binding Hamiltonian as well as some supplemental discussions on various related issues, including 
the band structure fitting, a comparison with previous theoretical studies, the one-loop RG equations, and the 
Luttinger parameters away from the condition of atomic orbital rotational symmetry. 


A. Solution of orbital quasi-degenerate tight-binding Hamiltonian and fitting the band structure 


The tight-binding Hamiltonian discussed in the main text is given by 


H 0 = - 




E fU) YMjO') i R2) Mm) j(m') ,(3) Mm) i(m') ,(4) Mm) ,(m') i 

t L mm ,Uj nIa a nl'cr ' L mm ,Uj nIa a nl'a ' L mm ,Uj nIa a nl'a ' L mm ,Uj nIcr n+1 /'a * 

E E + E E Et <si) 


z^ z^ 

n,a m,m' ///' 

■x 

n,cr m,m' I 


Where, d^l annihilates a Cr 3d-electron moving along the z-axis at the n-unit cell, with spin polarization a (=t,i), 
orbital component m (= 0 ,± 2 ), intra-triangle location /(= A, B,C), as well as the conjugate triangle location /(= 
A, B , C ). Here, the site locations correspond to (a, £) introduced in the main text as: A = ( 1 , 1 ), B = ( 2 , 1 ), C = (3,1), 
A = ( 1 , 2 ), B = ( 2 , 2 ), C = (3,2). The notation (a,£) introduced in the main text is explicit to accommodate the 
group representation , while (A, H, C) or (A, 5, C ) are more transparent. Either notations will be used for convenience 
in this SM. The hopping parameters t^ m , ( z = 1,2,3 and 4) are those for the n.n. sites shown in Fig A. In the 
present orbital quasi-degenerate model, these intraorbital and interorbital hopping parameters satisfy the relationship 
tmm' = Vfi™ fi\m\\m'\ for m 7 ^ m' with the ratio |? 7 | < 1 ( In the following we take 77 > 0 without losing the generality). 

= trn $ mm' is the next nearest neighbor (intraorbital) hopping along the tube direction. The CEF term is given 
by E m (or E m ), with n^ a (or n^_) being the corresponding density operators. 


In order to solve the non-interacting Hamiltonian for a single cluster [CrAs]6 in the n-th unit cell, Ho,m we introduce 


,(m) _ (Am) 
an ~ \^ a (n,z4)’ 


,(m)_ ,{m) 

a (n,A)’ in,BY 


d^ C y for each m = 0, ±1, ±2 (the spin index a is implied). Because of 

orbital mixing between m = ± 1 or ±2, we need to introduce another set of base [d^ ± d^R 171 ^]. Then 

all thirty atomic orbitals (MOs) in a unit cell can be described by dn \ with r — 0, ±1, ±2 ( denote dn ~= dl m_0 ^) 


The 


MOs of the n-th cluster, d ^ ( C ( n q,l)’ C (n!l,2)’ C \n,2,l)i ^ 7 ^ 2 , 2 )’ C \n,3 


S T ) 


(n,3,l) 5 (n,3,2) 


)• 


are defined as the base which 


diagonalizes Hamiltonian Ho, n = Cn^ H^Cn\ Here, H ^ = diag{E^ 7 ^A is a diagonal matrix with eigenvalues 


E\ 


The MOs can be obtained by Cn } = (R 0 Qo)dn } with uj = e Z( G ip = 27r/3, R and Qq are given in Eq.(3) 




7(0 


m 


7 


the main text. Explicitly, we have^ 


„( T ) 

-'KM) 

„( T ) 

'(n, 1,2) 

„( T ) 

"(n,2,l) 

(r) 

2 (n, 2,2) 

„( T ) 

"(n,3,l) 

„( T ) 


*(ra>A) ^ “(n.B) ^ (n,B) ^ u (n,C) 


_ J( T ) + JK _ J( T ) + J( r ) _ \ 

A) “(n,A) + a (n,B) a (n,B) + “(n,C) “(n.C) j ’ 


(n,C) 

t T ) 


~^J= (^(ji.A) + 4l4i + + <4l4l + 4lV'A + E 

— f J (T) 
ve l <"• 

+ < 4 A + wd (T) 

— ( J (T) 

V<"- 

— ( J (T) 

^ v< n - 


A) 


(n,A) 


(■ n,B ) + W ^(n!s) + W 1( ^(n!c) + W " 


A) “ <^A) + W ^(n!s) “ ^(nis) + W M^C) “ W ’ 


- 1 J (T) - 

“(n,C )) ’ 
“(n,C) y 1 


V 6 


} + d\ T J 3 , + LO 1 d { ^ B) + iO 1 d ( ^B) + ^(nlc) + W< ^(n!c)) > 


A) 


(n,A) 


\n,B) 


(n,C) 


(n,< 


C K3,2) = ^( d lniA)-^A)+ W ld '(n!s)- W ^B) + W< ^»!c) “ W ^C)) ' 


(5 2 ) 

(53) 

(54) 

(55) 

(5 6 ) 

(57) 


In order to solve the whole tight-binding Hamiltonian Ho, we need to introduce d^ = 
{$k\y tfk\y ^(/Tc)’ ^(aT^)) 5 an d the corres P on ding in the momentum space. Then, the 

Hamiltonian can be diagonalized by C^ r) = , in the form of 

Hq = S T k C'l^H q^C^. Here, H^ r k = diag{£^^(k)} is the diagonal matrix with eigenvalues £^^(k) given by 


£ (a,()( k ) = -S i!l) + e (I? 2 ) + (-^V^Si) “ e K 2 )) 2 + 4 ^ 


(S 8 ) 


In above, k is the crystal momentum along the tube direction, a£ t ^ = [1 + sign(r)\r]\]\ a for a = 1,2, 3, with Ai = 2 and 
A 2 = A 3 = — 1 the eigenvalues of H, e|^ = 4 -^ + 2tr cos ^~ ) Et p T = | 4 - 3 ^ + 4 - 4 ^e* fe |, and fi the chemical potential. 

Accordingly, the electron operators in MO bands are given by C^ = R® diag{Q^ a r \k)}d ( k \ Here, diag{Q^\k)} 
is a direct product of sub-matrices Q^ (fc) defined for each eigenstates A a of the C 3 rotation as given by 


Q ( a ] (k) = 


(r) 

,V> 

(t) jQ 
S1Ila (a,2) e 


cos a 


( T ) p-iOr 


r) 


(S9) 


A 4 ) sin k 


with tan 9 aT = (3) (4) 


(r) 

cos a; M = 


(a,2) / 

£< T > -e< T) 


— 1 an d A (r) = Q.o %A) 

iT+M'’ COS fc ' ( a ’£) ^/i_|_(a| T ^) 2 ’ ( a >£) Par 


(t) 

The Fourier transformation of C k ; back to the spatial space does not return exactly, though similar, to the forms as 
defined by Eq.(S2-S7), because the corresponding coefficients in each terms are now /c-dependent due to the intercell 
coupling along the tube direction. This feature will in general lead to various long-range electron correlations among 
the MOs of different unit cells. However, upon summation over the whole all unit cells, the slowly varying terms 
dominate the contributions. So as long as only the local interactions of MO bands are concerned, Qo or Eqs.(S2-S7) 
can be used to deduce these interactions as given in the next section. 

In order to fit the DFT band structure, we re-express the eigenvalues in a more explicit form 


£(a£) (fc) = COS k + /X + 


Et + En- 


[l+» 9 „(r),|A„ 1 1 ,„ +1?) + ( _ 1){ {t m _ { 


( 2 ) 


e t - e t 


[1 + sign(r)g\X ( 


-) 2 + 4|4 3) + t^e ik | 2 > .(S 10 ) 


As explained in the main text, we fit the three active DFT bands <a, /3, and 7 along the T-A direction ( the tube 
direction or the c-axis)^^. The non-degenerated 7 band is contributed mainly from the d z 2 orbital with m = 0, is 
fitted by the singlet MO band indexed by (r = 0 , a = 1 , £ = 1 ), with the energy s[°l(k), 


£( 7 ) (k) = — 2 £q 5 ) cos k + fi + 


En + Ed 


-4 1 ) -4 2) + \/(4 1 ) -4 2) - 


Ed — Ed 


) 2 +4|4 3 ) +4 4 ) e^| 2 . (Sll) 
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The a- and /3-bands, which are degenerate along the T-A direction, are fitted by the MO bands indexed by (r = 
—2, a = 2, £ = 2) and (r = —2, a = 3, £ = 2), with the energy £^ 2 ) (&) = ^( 32 ) W 


£(a,/3) (fc) = COS fc + /i + 


E 2 + E 2 


+ (! + ??) 


4 1} +4 2) 


1 


'(4 ) -4 ) + 


-E 2 — -Eh 

1 + Tj 


) 2 +4|4 3) +4 4) e ife |fS12) 


Both Eqs.(Sll,S 12 ) take the form as £(k) oc A T + B T cos k + ^/(7 T + D r cos 2 each with four independent coefficients 

T t , £? r , C T , and _D r . For the 7 -band, To = /x+ - E °+- e ° , E?o = -2tg 5 \ Cq = (44 ~4^ — ■ E °~^° ) 2 +4(4 3 ^—47 2 > 

Do = 164 !j f‘ 4 ’. The best fitting is given by To = —1.306 eV, £?o = 1-049 eV, Co = 2.194 eV, and Co = —2.121 eV. For 

[{t w_ t ^_E 2 ^ ) 2 +^-t^n 


the a- and /3-bands, A 2 = /x + 

x2-(3).(4) 


_ .. 1 E 2 +E 2 1 I +77 /A 1 ) 


+ i T 2 (4 i; +4 2) ) ) B 2 = -2t y °\ C 2 = 


— _o-/-( 5 ) C" — ( 1 + y ^) 


D 2 = 4(1 t\ • The best fitting is given by A 2 = —1.500 eV, H 2 = 0.962 eV, C 2 = 3.450 eV, and D 2 = —3.074 

eV. The negative D r implies opposite signs of t ^ and t^ • The relatively large value of C T implies a sizable difference 
in length or electron occupation between the two conjugated triangles as already indicated in the DFT calculations^ 
Of course, the above fitting is by no means rigorous, given the fact that the renormalization effect may be not 
adequately accounted in the DFT band structure. It is also possible to fit the DFT band structure within a reasonable 
approximation by other sets of parameters. For instance, we can use relatively smaller parameters F>( r ) and C( r ), but 
positive D r , the overall lineshape and band width are still closed to the DFT results. In comparison with the bare 
tight-binding parameters these fitting parameters should be all effective after renormalization. As far as the three 
active MO bands <a, /3 and 7 are concerned, the fitting formulae involve eight independent coefficients. Because the 
total number of free tight-binding parameters used in fitting exceeds eight, we cannot determine these parameters 
uniquely. On the other hand, the precise values of these parameters are not important in our present study. As we 
have shown in the main text, only the symmetry property of the active MOs and local interactions between them 
play the most crucial role in the resultant Luttinger liquid theory. 


B. Molecular orbital interaction parameters and comparison with previous theoretical studies 

In the main text of this paper, we have considered the microscopic atomic orbital (AO) Hubbard model with the 
local intraorbital Coulomb interaction V and interorbital Coulomb interaction [/', Hund’s coupling Jh , and pair 
hopping Jp, among all five atomic 3d-orbitals m = 0, ±1, ± 2 . Explicitly, the interaction matrices take the following 
forms 


/ U U' U' U' U' \ 

U' U U' U' U' 

U' U' U U' U' 

U' U' U' U U' 

V U' U' U' U' u ) 


/ 0 Jh Jh Jh Jh \ 
Jh 0 Jh Jh Jh 
Jh Jh 0 Jh Jh 
Jh Jh Jh 0 Jh 
V Jh Jh Jh Jh 0 / 


( 0 Jp Jp Jp Jp ^ 


J v 


0 Jp Jp Jp 


Jp Jp 0 Jp Jp 
Jp Jp Jp 0 Jp 

V Jp Jp Jp Jp 0 / 


(S13) 


Using the inverse mapping from the MOs to AOs, various two-particle interactions among the MOs are induced by 
the above local AO interactions. In general, the induced MO interactions are non-local due to the inter-cell hopping. 
As far as the local MO interactions are focused, the matrix Q 0 or Eqs.(S2-S7) can be used in deducing the local MO 
interactions given in Eq.(5) in the main text, as the short-range interactions are mainly due to the slowly-varying 
part. Such the local MO interactions depend on local AO interactions via various products taking the forms like 


AOi) AO 2 ) j(t 3 ) j(t 4 ) 

(n,a,£),cri (n,a,£),(T2 a (n,a,£),cr 3 a (n,a,£),cr4 


= £ 

i£,i 


[Qo 1 ®^ • [<5o 1 ® R 7x3,a 3 ,€ 3 ) • [Co 1 ® ^ 


C, 


4 On) 

(n,a 1 ,£1 ) 


c, 


rtto) 


t ti (n.flo., 


c: 


l( T 3) 




c 


r( T 4) 


(S14) 


In the previous section, by solving the tight-binding Hamiltonian and fitting the DFT band structure, the three 
active molecular orbital bands v = 1,2,3 are identified. Based on the symmetry argument, v = 1 corresponds 
to the 7 -band with m — 0, denoted by the molecular orbital electron annihilation operator ci(n) = cj^ ^ (with 
r = 0, a = 1, £ = 1) in Eq.(S 2 ). While v = 2 and v = 3 correspond to the a- and /3-bands with m = ±2 , denoted by 
the molecular orbital electron annihilation operator C 2 = c|~^ 2 ) r = — 2 ,a = 2 ,^ = 2 ) and C 3 = (with 
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r = — 2 , a = 3,£ = 2) in Eqs.(S5) and (S7), respectively. After some tedious but straightforward algebras thanks to 
symmetry properties of the matrices R and Q o, the following interaction matrices for various two-particle Coulomb 
interactions and Hund’s couplings for the above three active molecular orbitals are deduced: 


U 


J = 


Ux 

U 12 

Un \ 



U21 

u 2 

U23 = 


u' 

12 

U31 

U32 

U 3 J 


U' 

12 

0 

J12 

J13 \ 

( 

0 

J21 

0 

J23 = 


J H 

6 

■hi 

J32 

0 / 

\ 

J H 

6 


u 

6 


Jh 

24 

Jh 

24 


UL - Jh. 

12 24 

U-\-U' -\-Jh -\-Jp 
, 12 

U-\-U -\-Jh~\~Jp 


UL - Jh. 

12 24 

U-\-U +Jh+J p 
, 48 

U-\-U -\-Jh~\~Jp 


48 


12 


Jh 

6 

0 

U+U' + Jn+Jp 
12 


Jh 
, 6 

U-\-U -\-Jh~\~Jp 
12 

0 


In addition, a three-band interaction term 


1123 


c it( n ) c u( n ) c 2+( n ) c 3t( n ) + (2 3) + h.c. 


„t 


(S15) 


(S16) 


(S17) 


emerges, with J 123 = J p / 6 . Note that this peculiar interaction is absent in the previous theoretical studies*^. However, 
as we shall show in the next section, this term does not influence the instabilities we concern. 

Usually, we assume the orbital rotational symmetry of the interacting Hamiltonian in the atomic orbital Hubbard 
model, i.e., U' = U—2Jh and J p = Jh - 7 For a given material, it is understood that values of these local interactions can 
be meaningfully determined when band structure calculations as using the constraint density functional approaches are 
implemented by a complete description of the corresponding atomic Hubbard model. So far such detailed calculations 
on the compound K^CrsAss are not yet available. However, these local atomic interactions can be inferred from those 
in other Cr-based oxides, like SrCrC^. According to Ref4, one has U ~ 2.7 eV, Jh ~ 0.42 eV in SrCrCU. These values 
are relatively smaller than but still closed to those of Fe 2+ systems like SrFeCU or iron pnictides, where U ~ 3 — 5 
eV, Jh ~ 0.50 — 0.70 eV^ The ratio Jh/U of these systems does not change too much. Hence as a rough estimate, 
we assume the similar values of SrCrC^ for the present compound, obtaining the local molecular orbital interaction 
matrices from Eqs.(S15,S16): 

/ 0.45 0.14 0.14 \ / 0 0.07 0.07 \ 

U « 0.14 0.45 0.11 , 0.07 0 0.45 . (S18) 

\ 0.14 0.11 0.45 J \ 0.07 0.45 0 J 


Notice that it is always possible for these values to be fluctuated within 20% or even more for a given material by 
using different approaches. With this understanding, it is interesting to compare our results with those in Refs££. 
In our Q 1 D Luttinger theory, while the channel ”3” is always critical for a much wider regime U > Jh , the channel 
”2” is expected to be located at the spin triplet phase because Jh/U ~ 0.16. However, this triplet state should close 
to the border of the SDW phase, Jh/U = 0 . 2 , as shown in Fig.3 in the main text. 

In Ref£, the following interaction matrices for phenomenological molecular bands are introduced: 


while in Ref£, 


(u x U 2 U'\ 

Uz = \ U 2 U x U^ , J z 

V ^2 ul u{ J 


0 j j'\ 

J 0 J' , J p = J z , 
J' J' 0 


(U V v\ 

u w = [ v U V , Jw 
\V V u J 


0 Jh Jh\ 

Jh 0 Jh , J P = Jw- 
Jh Jh 0 J 


(S19) 


(S 20 ) 


Here, we use the same notations introduced in Refs££, respectively. The relationships U[ — Ui, U' 2 = U 2 , and 
0.5 < J'/ J < 2.0 were used in calcuations^. It should be understood that (i) the parameters used in these studies are 
for the molecular orbital bands; (ii) they do not rigorously correspond to our interaction matrices given by Eqs.(S15- 
S17). Nevertheless, we may expect that if the largest elements in the interaction matrices U and J dominate in the 
random phase approximation (RPA) approach, the values of atomic interactions estimated in our case correspond to 
the case of U\ = 0.45, J = 0.45, and J/U\ = 1 in Ref£. So the channel ” 2 ” should correspond to the regime with 
small U but large J/U in Ref£. As shown in Fig. 3(a) in Ref£, this phase favors the spin triplet fy^^-y 2 ) pairing 
state. Notice that these values are out of the calculated regime of the phase diagram Fig ( 6 ) in Ref^, but the same 
tendency is reached if this phase diagram is extended to large ration of J/U. 
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C. One-loop RG equations and absence of the three-band instability 


Owing to the symmetry between the MO bands v = 2 and v = 3, there are twelve permissible perturbations which 


can be classified in terms of the formal g-ology descriptio n 10 ! 11 using the same notation : 


01-14. 


•Hint = 1 -dOw 

<J<y' 

+ g^J2 R l R ^ L l'Lia' 

aa' 

9^ ^L^2a'^3 a '^3a + ( 2 ^ 3) 

aa' 

- 5 i 2) ^i?LLv4^3, + (2o3) 

aa' 

+ #2 ^ ^ 2 ( 7 ^gLzq'Rsa' + ( 2 ^ 3 ) 

aa' 

+ #2^ ^ ^3a / +(2^3) 

aa' 

+ #4^ + (2 3) 

aa' 

+ #4^ y^i?2 C r^2cr^2cr / ^ / 2cr / +(2^3) 

aa' 

- /« ^[4 CTJ R la ,0,X 2a + (2 ^ 3)] + (iZ L) 

aa' 

+ /< 2 > + (2 ^ 3)] +(R^L) 

aa' 

+ + (2 3)] + h.c. 

aa' 

+ ^ y^[-^io-^i 0 -/-^2o- / ^3g + (2ty 3)] + h.c.. (S21) 

aa' 


In this notation, the degeneracy between the MO bands v = 2 and v = 3 are explicit. The one-loop renormalization 
group (RG) equations for these coupling constants evolving with increasing scaling parameter l are of common type, 
given 


dig {1) 

dig {2) 

dig { i ] 

dig ( i ] 

dig { 2 ] 

digi 2) 

dig { A ] 

dig { 4 } 

dif {1) 

dif {2) 


—2(<7 (1) ) 2 - 4 uv 
-(<? (1) ) 2 - 2 u 2 - 2v 2 

—2(<0) 2 ‘-g ( i i g ( 2 ) -2mu 

- 25 ^ 1) ^ 1) - g {2) g {2) + <7i 2) i/i 2) -u 2 -v 2 

-25i 1 ) 5i 2) - 2<?2 1) fl2 2) + 5 i 2) 54 1} “ 4fl 2 1) <?4 1) + 2 <? 2 1) 54 2) “ 2 «£ 

-Cffi 15 ) 2 - (ffPb 2 - c ^ x) ) 2 - a 2 - v 2 
2^ 2) 0 ) -2(0 ) ) 2 -2(0 ) ) 2 
(#) 2 - G^) 2 

—2(/( 1 )) 2 + 2uv — 2v 2 

-(/ (1) ) 2 +u 2 

( 2 /! 2 ) — < 7 ^ — ^! 2 ! — g^)u — (/ (1) + flb + 9^)v 

-(-2/« + 5 (1) + sb + #)« - (4/ (1) - 2/( 2 ) + ^ + 5 (2) + ~di 2) )v • 
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U/27tv f U/2ttv f 


FIG. SI: The dominating perturbations in the one-loop RG equations. (a)All the initial parameters are given by Eqs.(S15- 
S17,S22); (b) The initial values of u and v have the opposite sign using v/2itvf = Jh/6 — 0.02 and u/2ttvf = Jh/6 for small 
Jh, while the other parameters are still given by Eqs. (S15-S17,S22). 


In above, g = 2 ^ vf , and 


5 (1) = Ui, 

~(1) _ T 

91 — ^23, 

9 { 2 ] = 2C/23 + y 23 , 

g { J = u 2 , 

/ (1) = J 12 , 

^ = J\23i 

The corresponding bare MO interaction parameters, which are related to initial values of perturbations defined in 
Eq.(S22) in solving the RG equations, are given in Eqs.(S15-S17). The initial values of perturbations u and v are the 
same in the present model, both induced by the three-band interaction J 123 . When l approaches a sufficient large 
cutoff scale /*, the generic asymptotic solutions take the form gf^ oc (/* — An instability takes place whenever 

some of the couplings diverge at a finite scale length /*. We numerically determine the most divergent coupling when 
l approaches /*, which is taken to be unit, from below. For usual 3d-transition metals, 0 < Jh/U < 1, we find that 
gi 1 ) or g^ dominates in the large or small Jh/U regimes, respectively, as plotted in Figure ISTt ai. These two regimes 
correspond to the known single-band and the two-band instabilities respect ivelyi^r— . 

It was suggested that in addition to these one-band and two-band instabilities, there may be a new instability 
driven by the presence of all three bands due to the perturbation u or v. We have checked that this instability does 
not occur if the initial u and v have the same sign as determined in the present case. We have also checked that 
if the initial values of u and v have the opposite sign, say, assuming v/2i tvf = J#/6 — 0.02 and u/2ttvf = Jh/ 6 
for very small Jh , there is a regime where v dominates over all other perturbations as shown in Fig. ISTl b). It 
is interesting to recall that the electron-phonon coupling, which has not been adequately considered in the present 
study, may result in deviations of effective perturbations away from the initial values determined in Eq.(S22). How 
such modification upon electron-phonon coupling takes place and whether it influences the superconductivity in the 
K 2 CrsAs 3 compound deserve further investigations. 


g {2) = Ui, 

9 ^ = J 23 , 

gi 2) = 2U 23 + | J 23 , 

gf ] = U 2 , 

/( 2 ) = 2 u 12 + | J12, 

V = Jl23- 


(S22) 


D. The Luttinger parameters away from the condition of atomic orbital rotational symmetry 

The rotational symmetry in the interacting part of a given system is frequently assumed in literatures, by using 
the relationships J p = Jh and U' — U — Jh — Jp at the level of atomic orbitals^. The relationships are adopted 
in plotting the phase diagram Fig.3 in the main text. It is possible that in realistic and complicated systems this 
symmetry could be broken, or the above relationships could not be respected. In order to understand whether our 
results are still valid in this case, we consider a small deviation AU away from the rotational symmetry, by assuming 
U' = U + AU — 2 Jh- It is straightforward to show that the diagonal channels are similar to the case with AU = 0 as 
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given by Eq.(7) in the main text, while the model parameters appear in the Tomonaga-Luttinger Hamiltonian Eq.(8) 
are given by 


V F 

. Ui 

tc2 — 

V F 

. u 2 

2Ui2 


2C/23 

2 + 2tt’ 

2 + 2tt’ 

— 5 

7 T 

be — 

7T 

v F 

Ui 

t s 2 = 

v F 

U 2 

J 12 


J23 

2 2t r’ 

2 2n ’ 

a a =-, 

2tt ’ 

b s = 

2i r 


Then, the Luttinger parameters are given by K 1 i = J 2 \ ? with 


A 7 ,i — -(tj i + t 7 2 ) + - + yj8a ^ + (b 7 + t 1 2 — ^ 7 i) 2 ^ , 

A 7j 2 = ~(^ 7 i T ^ 2 ) T 2 ^7 + (6 7 T £ 7 2 Lyi)^ 5 

A 7j 3 — t 7 2 bry. 


Explicitly, we have three pairs of charge and spin Luttinger parameters in the corresponding channels: 


Where, 


K Ct 1 = 

= 


1 

\/l + Gi + G2 

1 

Vi — Gi + G3 


i^c,2 = 

K S ,2 = 


1 

a/1 T Gi — G 2 

1 

Vi — Gi — G3 


^c,3 = 1, 

^ S ,3 = I- 


G\ 


G 2 


g 3 


u 


+ 


At/ 


4^^ 127TVF 



'U + AU _ 5 J H V 

127™f / 



\ 12^^ 


C/ + AEA 2 

127Tt>_F / ’ 


(S23) 


(S24) 


The fact that the channel-” 3” is still in the critical phase with K s?3 = K c ^ = 1 is apparently due to the degeneracy 
of the MO bands v = 2 and v — 3. Based on these expressions, we find that all our results in the main text remain 
unchanged because the role of A U is to modify the value of U in a simple manner. 


Electronic address: fxyong@163.com 
Electronic address: daijh@zju.edu.cn 

1 H. Zhong, X.Y. Feng, H. Chen, and J. Dai, Phys. Rev. Lett. 115, 227001 (2015). 

2 There could be a prefactor (—1) for the sectors A a = — 1 (for a — 2, 3) related to the sector A a = 2 (for a — 1) as defined in 

Eq.(S9). This global prefactor can be absorded into the base definition and does not infleucne the final results. 

3 H. Jiang, G. Cao, and C. Cao, Sci. Rep. 5, 16054 (2015). 

4 X. Wu et al., Chin. Phys. Lett. 32, 057401 (2015). 

5 Y. Zhou, C. Cao, and F.C. Zhang, arXiv: 1502.03928. 

6 X. Wu, F. Yang, C. Le, H. Fan, and J. Hu, Phys. Rev. B 92, 104511 (2015). 

7 C. Castellani, C. Natoli, and J. Ranninger, Phys. Rev. B 18, 5001 (1978). 

8 L. Vaugier, H. Jiang, and S. Biermann, Phys. Rev. B 86, 165105 (2012). 

9 H.J. Xiang, S.-H. Mei, and M.-H. Whangbo, Phys. Rev. Lett. 100, 167207 (2008). 

10 J. Solyom, Adv. Phys. 28, 201 (1979). 

11 T. Giamarchi, Quantum Physics in One Dimension, Oxford University Press, Oxford, UK, 2003. 

12 Y.A. Krotov, D.-H. Lee, and S.G. Louie, Phys. Rev. Lett. 78, 4245 (1997). 

13 J. Gonzalez and E. Perfetto, Phys. Rev. B 72, 205406 (2005). 

14 D. Carpentier and E. Orignac, Phys. Rev. B 74, 085409 (2006). 


























